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Abstract. Noise-assisted transport in quantum systems occurs when quantum 
time-evolution and decoherence conspire to produce a transport efficiency that 
is higher than what would be seen in either the purely quantum or purely 
classical cases. In disordered systems, it has been understood as the suppression 
of coherent quantum localisation through noise, which brings detuned quantum 
levels into resonance and thus facilitates transport. We report several new 
mechanisms of environment-assisted transport in ordered systems, in which there 
is no localisation to overcome and where one would naively expect that coherent 
transport is the fastest possible. Although we are particularly motivated by the 
need to understand excitonic energy transfer in photosynthetic light-harvesting 
complexes, our model is general — transport in a tight-binding system with 
dephasing, a source, and a trap — and can be expected to have wider application. 
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Recent experimental studies of photosynthetic light-harvesting complexes have 
confronted us with the fact that at least some of these systems exhibit excitonic 
coherence that is surprisingly long considering their noisy environment [TJ [2j [3] . This 
makes it clear that if we are to understand their high light-harvesting efficiency, we 
must study the ways in which quantum transport is affected by the interplay of 
coherence and noise [1 G2 EJ 1 1 HS CCD HI HI HI H3 H3 HI HE1 Q33 \M HQ ■ It has 
been found that noise can enhance quantum transport in model excitonic Hamiltonians 
El [5] , a phenomenon called environment-assisted quantum transport (ENAQT) or 
decoherence- assisted transport. 

In the simplest approach, different environments around each chromophore lead 
to a tight-binding model with sites that have different energies (disorder). Because 
of disorder, the exciton becomes localised through coherent phenomena such as 
destructive interference or Anderson localisation [5J [3 [51 HI] . ENAQT is then simple to 
understand: noise can destroy the coherent localisation, helping the exciton reach the 
trap site and increasing the efficiency. Alternatively, decoherence has been described 
as fluctuations of site energies which can transiently bring levels into resonance, 
facilitating transport. 

If these interpretations were the whole story, ENAQT would be impossible in 
ordered systems, those without energetic disorder. The absence of ENAQT in ordered 
linear chains was predicted at least twice: for example, Cao and Silbey predict "the 
lack of environment-assistance in linear-chain systems" [T2] while Plenio and Huelga 
note "the expectation that noise does not enhance the transport of excitations" in 
uniform chains [5 . This expectation is strengthened by proofs of the impossibility 
of ENAQT in end-to-end transport in ordered chains HI]- We revisit the ordered 
chain and show that the case of end-to-end transport is the only case where ENAQT 
is impossible: its absence in end-to-end transport is the rare exception that we are 
able to explain. We anticipate that these findings will shed light on the efficiency 
of transport in ordered excitonic systems, whether artificial, such as J-aggregates, or 
natural, such as the LHII complex in purple bacteria and the chlorosome in green 
sulphur bacteria. 

Results related to ours were reported by Gaab and Bardeen [6 J . They considered 
ordered systems and noted that the environment can sometimes enhance the "effective 
trapping rate." By contrast, we focus on the trapping efficiency, a measure of how 
often the exciton is productively trapped as opposed to lost, regardless of how fast 
the transport is (loss is not modelled in [6], so the efficiency is always I). We do 
not average over initial sites, which allows us to explain the absence of ENAQT in 
end-to-end transport. 

1. The model and the definition of ENAQT 

The ordered system we consider is a one-dimensional array of N identical sites, coupled 
to their nearest neighbours and described by the Hamiltonian 

N-l 

H = VJ2 (\m)(m + l\ + \m + l)(m\), (I) 

where V is the coupling strength and h — 1. This Hamiltonian is equivalent to a 
system of coupled two-level systems that is restricted to the single-excitation sector. 
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Figure 1. Definition of ENAQT. a) Three-site chain with initial site |2) and 
trap site b) Transport efficiency r) in the three-site chain, with trapping rate 
k = O.f and loss rate /i = 0.01. r\ is maximised at the optimal dephasing rate 
7opt = 0.319. ENAQT is the magnitude of the enhancement, £ = ??max — To, or > 
in this case, § = 0.038. 



To study the efficiency of transport mediated by this Hamiltonian, we introduce 
two distinct attenuation mechanisms. First, the particle is irreversibly lost from each 
site at an equal rate /i, modelling processes such as exciton recombination. Second, at a 
particular trap site |r), the particle can be trapped at a rate k, modelling, for example, 
the transfer of an exciton to a photosynthetic reaction centre. These attenuation 
mechanisms are incorporated by adding a non-hermitian part to the Hamiltonian, 

-fatten = \ m ) (H ~ iK \ T ) ( T \ ■ ( 2 ) 

rn 

Loss and trapping both result in particle disappearance and have the same 
mathematical form; the distinction is that we consider the energy carried by lost 
particles to be unavailable while the trapped energy to be productively useable. The 
norm of the state at time t is the probability that the particle will survive that long. 

The attenuation mechanisms continuously reduces the particle's survival 
probability, so that after a sufficiently long time, t /J -1 , the probability of finding 
the particle is negligible. If p{t) is the system's density matrix at time t, the probability 
of trapping the particle in the interval [t, t + dt] is 2k (t \p(t) \ r) dt. The efficiency of 
transport, for initial state /?(0), is then the overall trapping probability, 



/■OO 

t] = 2k (r\p(t)\r) dt. 
Jo 



(3) 
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Likewise, the probability of loss is rf = 2/i^ m (m\p(t)\m) dt = 2p f^°trp(t)dt 
and these branching ratios satisfy rj + rf = 1. 

Environmental effects are modelled as (Markovian) pure dephasing, acting 
independently on all sites with an equal rate 7. We choose dephasing because it 
is one of the simplest forms of noise, giving us a single-parameter minimal model 
that exhibits the desired behaviour. Insofar as dephasing is the appropriate limit of 
several more realistic noise models, we can expect qualitatively similar effects if more 
complicated environments are considered. The dephasing superoperator T> is defined 
through {'Dp) nm = ^ 27 (1 — 5 nm ) p nm and the resulting complete equation of motion 
is 



where the total Hamiltonian is H = Hq + H attcn . We note that this master equation 
has been solved exactly for the case k = p = both in the single-particle situation 
(55J 123] and in the non-equilibrium setting for many particles [24, 25j. In the following, 
we give a method for the exact (single-particle) solution for any chain length and any 
k and p, allowing us to calculate the efficiency. 

The coupling V sets the energy scale, so we can take V — 1. Then, at every choice 
of loss and trapping rates p and n, the efficiency r\ is a function of the dephasing rate 7, 
see Fig. [TJ). We observe that if 7 is very large, the particle will be localised at its initial 
site due to the Zeno effect. Therefore, it will not be able to reach the trap before it is 
lost, meaning that 77 — > as 7 —¥ 00. Consequently, the maximum transport efficiency 
?7max will occur at a finite 7 op t > 0. ENAQT occurs if 7 opt > 0, and is defined to be 



We can now consider efficiency and ENAQT as a function of k and p. The 
following descriptions are all borne out in the example in Fig. [51 which shows 770 , £, 
and 7o P t as a function of k and p in the finite system of three sites with the trap at one 
end and the initial site in the middle (see Fig. [TJi). Regardless of the number of sites, 
several limiting cases can be easily understood. First, if k <C p. the particle will be 
lost before it can be trapped, regardless of the amount of dephasing present. Second, 
at large k, the Hamiltonian term — in |r) (r| presents a high potential barrier for the 
particle, meaning that it will be largely unable to access the trapping site (this is the 
case even though the potential is imaginary). Consequently, high trapping efficiency 
is possible only in the regime of small p and intermediate k (see Fig. [5Ji). High 
ENAQT is also only possible in this region because outside of it, loss is so dominant 
that dephasing will be unable to appreciably increase trapping. This is illustrated in 
Fig. dbi where it can be seen that ENAQT is large only where 770 is neither very close 
to nor very close to I. That is, ENAQT occurs when neither trapping nor loss is 
severely dominant, meaning that noise can push the balance in favour of trapping. 

2. ENAQT on a finite chain 

2.1. Analytical solution 

Although there is no general solution, ENAQT can be determined analytically in every 
particular finite system, which is how Plenio and Huelga proved that £ = in the case 
with the origin and trap at opposite ends of the chain [5], The solution is by Gaussian 




(4) 



£(Mj«) = l)max - V0- 



(5) 
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Figure 2. Characterisation of ENAQT in the three-site system from Fig. 
Each of the three parameters defined in Fig. [Tp (ryo, £, and 7opt) is displayed as 
a function of trapping and loss rates, a) High efficiency in the regime with no 
dephasing is possible only for small loss and intermediate trapping. The two lines 
indicate the region where the efficiency is around 1/2, meaning that neither loss 
nor trapping dominates. They are computed from Eq. [6]('?o = a o/Po)- b) Large 
ENAQT occurs in regions not far from the lines in a, where dephasing can push 
the balance between loss and trapping in favour of trapping, c) In the lower left, 
( max occurs when 7opt ~ 0.4151. In the upper left, strong trapping presents a 
high potential barrier, meaning that strong dephasing is necessary for ENAQT. 



elimination (see Appendix), meaning that r\ is a rational function of 7, k, and /i. For 
the three-site example in Fig. [T^,, 

_ a 2 7 2 + Q07 + a o / fi \ 
V ~ fol 3 + + 0i7 + A) ' U 

where a?2 = 4k/i, ct\ = n (2 + 2n[i + 8/x 2 ) , — k (2k/i 2 + k + 4/r 3 + 2/i) , 

f3 3 = 8h 2 (k + fj,), /3 2 = 4/i(2 K 2 /i + K (8/ 1 2 + 3) +6/x 3 +4/i), ft = 2/sV + 
2k 2 pi (9/i 2 + 5) + 2k (20/i 4 + 20a* 2 + l) + 6 (4/i 5 + 6/r 5 + ,u) , and /3 = 2k 3 (/i 3 + fi) + 
k 2 (10// + 13a* 2 + l) + K/i (I6(i 4 + 29fi 2 + 6) + 4/i 2 (2/i 4 + 5m 2 + 2) . ENAQT is 
calculated by maximising this function with respect to 7. In particular, it can be found 
that drj/dj can only equal zero if -K 4 /i+4K 3 /i 4 -2K 3 /r 2 + 2K 3 + K 2 (20/i 4 + 7fi 2 + 9) /i+ 
32k/i 6 + 16k/i 4 + 7k/i 2 — k + 16/i 7 + 8/i 5 — 4/i 3 — 2/i > 0, meaning that there is a region 
in the (ju, k) plane in which ENAQT is impossible, as shown in Fig. In all other 
cases, ENAQT is strictly positive. Maximum ENAQT is £ max = 7 - 4^ « 0.0718, 
obtained as k and /i simultaneously tend to zero while keeping /i = k/2s/3. In that 

limit, 7opt tends to yj a 0.4151. 

"2.^. Limit of small attenuation 

It is difficult to form a simple, intuitive picture of ENAQT in this system that remains 
valid in all parameter regimes, and particularly when time scales converge, e.g. at 
K ~ 1 or fj, ~ 1. Nevertheless, there is a simple expression for ENAQT in the limit 
k, fj, <^ 1. In that case, both attenuation mechanisms are weak and can be treated as 
perturbations on the dephased quantum dynamics. In particular, because attenuation 
is slow compared to the quantum dynamics, we assume that we can only consider the 
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Figure 3. ENAQT in the three-site system from Fig. in the limit of small 
attenuation, with k = = 10~ 3 . As 7 — >• 0, the time evolution is merely perturbed 
coherent evolution, and the efficiency approaches the predicted value of i (dash- 
dotted line). As 7 — > 00, the Zeno effect ensures r] — > 0. In the intermediate 
regime, 7 ~ 1, the dephasing creates a fully mixed state, for which the predicted 
efficiency is j (dashed line), giving an ENAQT of ^j. 



average site populations in calculating loss and trapping. In the case with appreciable 
dephasing, the state of the system will quickly reach a completely mixed state, meaning 
that each site will host 1/N of the remaining population. In particular, the rate at 
which the particle will be trapped at the trap site will equal k/N. Similarly, all 
population will be lost at a rate /1, giving the efficiency 

_ k/N _ 1 m 

^ephascd ~ K / N + fi - 1 + Nfl/K V) 



In the coherent case, the system eigenstates are = y^ jv+i sm n+i w ith eigenvalues 
Afc = 2 cos j^q-fj-. Therefore, the amplitude of site I given an initial site to is 

Ulm{t) = 2 A ^ e -, t , cos ^ (g) 

u JV + l^ JV+1 iV+1 w 

3=1 

which can be used to show that the average population P; m = limy^oo h f Q \Ui m (t) \ 2 dt 
equals 

1 ( 1 „ 1 



P lm = N+1 ( 1 + 2 Slm + 2 S l,N+l-m J , (9) 

where Si m is the Kronecker delta function. From there we have ^coherent ~ 
(l + P^/x/^) • For ENAQT, we do not consider the case I — to, meaning that 
there are two situations, depending on whether the trap is opposite the initial site. 
Because partial recurrences can refocus excitation from the initial site to the opposite 
site, if the opposite site is the target, I = N + 1 — to, the average target population 
exceeds -4 and ENAQT is impossible. This explains, at least in the limit of small 
k and fi, Plenio and Huelga's observation of the absence of ENAQT in end-to-end 
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transfer. In transport between sites that are not opposite each other, 



1 + Nh/k 1 + (N + 1)h/k ' v ' 

Notably, £ depends only on the ratio (j,/k. The validity of these approximations is 
demonstrated in Fig. [3] The expression is more accurate for large ju/k because loss, 
by lowering all the amplitudes simultaneously, perturbs the time evolution less than 
trapping, which affects only one site. 



2.3. Other patterns 

Several patterns emerge in longer chains and when the locations of the trapping site 
and the initial site are varied. Table[T]shows the maximum possible ENAQT in chains 
up to N — 8 with all possible combinations of initial and trap sites. Each entry 
is calculated by analytically solving the equations of motion and maximising £ as a 
function of k and [i, as discussed above for the N = 3 chain. As we proved above, we 
can see that ENAQT is possible in all configurations except when the initial site is 
located opposite the trap. Furthermore, maximum ENAQT increases with increasing 
N. This is the opposite of the trend predicted by Eq. [TUJ and occurs because the 
high values observed in the table generally occur in the regime of small h/k, where 
the estimate of ^coherent fails. 

It can also be seen that, regardless of the trap site, £ max is equal for situations with 
initial sites m and N + 1 — to. This is the case even though ^(fi, k) is not equal in the 
two situations in general. In the limit of infinitesimally small k, /i, and 7, where the 
equality obtains, coherent time evolution proceeds before any appreciable dephasing 
or loss takes places. Therefore, coherent recurrences can occur, and although the 
recurrences are imperfect because the eigenvalues are incommensurable, they 
become arbitrarily close to perfect after a sufficiently long time. In particular, a 
particle initialised at to will refocus (arbitrarily close to perfectly) at N + 1 — to after 
sufficient time. On the longer time scale of loss, the two initial conditions become 
indistinguishable . 



2.4- High ENAQT in symmetric situations 

As shown in Table [TJ a high ENAQT of £ max = h occurs if the trap is in the middle 
of a chain with an odd number of sites, regardless of the initial site. In that case, 
the full Hamiltonian H commutes with the inversion operator P, defined as P \ j) = 
\N + 1 — j). The initial site |r) can be written as an equal superposition of symmetric 
and antisymmetric states, |r) = + \A)), where|S) = ^ (|r) + \N + 1 - r)) 

and \A) = -^(\t) - \N +1-t)). Because \A) is odd, P\A) = -\A), and H 
commutes with P, \A) remains odd under time evolution. Since the trap site is in 
the middle, it is even, meaning that the \A) component of |r) never gets mapped to 
the trap site and can therefore not be trapped. By contrast, the \S) component does 
get trapped. Therefore, the efficiency at zero dephasing is rjo = |. When dephasing 
is non-zero, the phase coherence in \A) is lost, meaning that now the particle can be 
completely trapped. In particular, if k S> (J., a negligible amount will be lost, meaning 
that r/ max = I, giving £ = r/ max - r/ = \. 
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Table 1. Analytically calculated maximum ENAQT in a chain with N sites, 
depending on the initial and trap (t) sites. ENAQT is not defined if the initial 
and trap sites coincide (X). All the values with the same N that are equal to three 
decimal places are exactly equal. 0.072 = 7 — 4\/3 and the rest are roots of more 
complicated, but known, polynomials. 



3. ENAQT on a circle 

If the transport takes place on a circle instead of a finite chain, the behaviour is 
qualitatively the same. We consider here the same situation as above, except that in 
Eq. [TJ there is an additional term coupling sites |1) and \N). 

In the regime of weak attenuation, k, fi <C 1, as in the chain, ?7dcphascd ~ 
(1 + N(x/k) 1 . In the coherent case, the eigenstate amplitudes are Uj = -l= e 27Tl 3 k / N 
with eigenvalues = 2 cos , from where the average population is 

- = (^(N(l + 5 lm )-l), iVodd 
lm \^(N(l + 5 lm + 8 Lm+N/2 )-2), iVeven. 

Because P; TO > -k if I = m + ^, there is no ENAQT if the initial site lies on the 
opposite side of the circle to the trap. In other cases, 

^~ 1 + Nh/k ~ 1 + N 2 ^/(N - 1)k > °" (12) 
As in the chain, this expression is more accurate for large /x/k. 
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Figure 4. The infinite system (top) and its characterisation (bottom). As in 
the finite chain (Fig. [2}, high rjo occurs for small loss and intermediate trapping, 
ENAQT is highest at intermediate 770, and the dephasing required for ENAQT 
grows rapidly with increasing k. 



Regardless of the number of sites, maximum possible ENAQT on the circle is 
JO, initial and trap sites opposite 

(max = i 1 . (13) 

I 2, otherwise. 

The situation is much simpler than in the chain, where Table [1] is needed. The high 
value of £ max occurs because the circle always has the inversion symmetry required in 
Sec. I 



4. ENAQT on an infinite chain 

ENAQT also occurs — albeit by a different mechanism — in infinite ordered systems. 
The dephased dynamics of a particle on an infinite chain is well understood in the 
absence of trapping and loss |261 1271 [T7| . In the fully coherent case, 7 = 0, the 
dynamics is ballistic, while increasing noise converts it to diffusion. 

To observe ENAQT, we introduce loss /i everywhere and trapping k, on all sites 
to the left of the initial site (see Fig. m,). In the coherent case, there is a sizeable 
probability that the particle, owing to its ballistic motion, will move far to the right, 
eventually being lost. In the opposite extreme, 7 = 00, the particle is completely 
localised at the initial site by the Zeno effect, and therefore eventually lost. In 
the intermediate region, ENAQT is possible because decoherence slows down the 
spreading sufficiently to prevent the rightward-moving component from escaping, but 
not strongly enough to prevent the particle from diffusing into the trap region. 
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Fig. [Dd shows numerically computed ENAQT on the infinite chain. The infinite 
line was represented by sufficiently many sites A to avoid the particle reaching the 
edges. Because the simulation time and A both scale as p~ l , the computation becomes 
expensive for small p, and we have imposed a lower cutoff of p = 0.1. The largest 
ENAQT found is = 0.1, k = 6.3) = 0.1233. Similar results are obtained with 
initial sites further from the trapping region. In those cases, ENAQT is smaller because 
the particle has to travel farther to the trapping region, but it remains finite for small 
p and intermediate k. ENAQT tend to zero as p — > 0, as p — > oo, and as k — ¥ oo for 
the same reasons as in the finite chain. A question we are not able to answer with 
numerical simulations is whether ENAQT is ever exactly zero or merely approaches 
zero asymptotically in the appropriate limits. 

5. Conclusions 

We have shown at least two different mechanisms for ENAQT in an ordered system. 
In the finite lattice with small k and p, it is caused by the fact that a dephasing- 
induced mixed state is more likely to be found at the trap site than a coherently 
propagated initial state, except if the initial and trap sites are opposite each other. In 
the infinite lattice, it occurs when dephasing slows down otherwise ballistic transport 
and prevents a portion of the particle from escaping far from the trap region. We 
leave open the questions of whether the various mechanisms of ENAQT (including 
the ones in disordered systems) can be understood in a unified picture and how they 
are influenced by different kinds of noise, of which pure dephasing is only one limit. 
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Appendix: Analytical calculation of the efficiency 

The initial state of the system p(0) can be understood as a vector in a Liouville space 
of dimension iV 2 . In order to calculate the efficiency, we aug ment it to p(O) with N 2 + 1 
entries. The first A*" 2 entries we call the state sector, and they equal p(0), while the final 
entry, the accumulator, we initialised to /5jv 2 +i(0) = 0. The Liouvillian C is likewise 
modified to C, an (A 2 + l) x (A 2 + l) matrix, where the top-left A 2 x A 2 elements 
equal £, and the remainder are set to 0, except for the entry £jv 2 +i,t' = 2k, where 
t' = 1 + (t — 1)(A + 1) is the coordinate of the population of the trap site r. That is, 
C couples the population of the trap site to the accumulator (but not vice- versa) with 
strength 2k. Because C does not couple from the accumulator to the state sector of 
pit), the time evolution of the state sector under C is equal to the time evolution of p(t) 
under C. During the evolution, the accumulator increases precisely at the rate 2np T i(t), 
meaning that /5/v 2 +i(oo) = ry, while the remaining elements of p(oo) are all reduced to 
0. In principle, one could calculate r\ by calculating p(oo) — lim^oo e ct p(0), but this 
appears to us to be too difficult analytically. 
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Instead of solving the initial-value problem, we solve a related steady-state 
equation. We begin by modifying £ to £ £ , which is the same except for C £ N2+1 jv 2 +i = 
e. Then we solve, analytically by Gaussian elimination or otherwise, the linear system 
of equations 

t E ~a = ep(0). (Al) 

Here, a is the steady state in the situation where p(0) is being injected into the system 
at rate e. In particular, since total probability being injected into the system is e, the 
fraction erj must go to the accumulator. Since the accumulator accumulates at a rate 
2k5v') we must have 5y/ = ej]/2k. Now, the accumulator component of Eq. IA1I is 
2na T > + etrjv 2 +i — 0, from where we can conclude that ct W 2 +1 = 77, meaning that the 
efficiency can be read out of the solution a . The result is independent of e, meaning 
that the procedure remains valid in the limit of infmitesimally small e, where Eq. IA1I 
reduces to a true steady-state equation, a = 

References 

[1] Gregory S Engel, Tessa R Calhoun, Elizabeth L Read, Tae-Kyu Ahn, Tomas Mancal, Yuan- 
Chung Cheng, Robert E Blankenship, and Graham R Fleming. Evidence for wavelike energy 

transfer through quantum coherence in photosynthetic systems. Nature, 446:782, 2007. 
[2] G Panitchayangkoon, D Hayes, K A Fransted, J R Caram, E Harel, J Wen, R E Blankenship, 

and G S Engel. Long-lived quantum coherence in photosynthetic complexes at physiological 

temperature. Proc. Natl. Acad. Set., 107:12766, 2010. 
[3] Elisabetta Collini, Cathy Y Wong, Krystyna E Wilk, Paul M G Curmi, Paul Brumer, and 

Gregory D Scholes. Coherently wired light-harvesting in photosynthetic marine algae at 

ambient temperature. Nature, 463:644, 2010. 
[4] Patrick Rebentrost, Masoud Mohseni, Ivan Kassal, Seth Lloyd, and Alan Aspuru-Guzik. 

Environment-assisted quantum transport. New J. Phys., 11:033003, 2009. 
[5] M Plenio and S Huelga. Dephasing-assisted transport: quantum networks and biomolecules. 

New J. Phys., 10:113019, 2008. 
[6] Kevin M Gaab and Christopher J Bardeen. The effects of connectivity, coherence, and trapping 

on energy transfer in simple light-harvesting systems studied using the Haken-Strobl model 

with diagonal disorder. J. Chem. Phys., 121:7813, 2004. 
[7] P Rebentrost, M Mohseni, and A Aspuru-Guzik. Role of quantum coherence and environmental 

fluctuations in chromophoric energy transport. J. Phys. Chem. B, 113:9942, 2009. 
[8] F Caruso, A W Chin, A Datta, S F Huelga, and M B Plenio. Highly efficient energy excitation 

transfer in light-harvesting complexes: The fundamental role of noise-assisted transport. J. 

Chem. Phys., 131:105106, 2009. 
[9] Alexandra Olaya-Castro, Chiu Lee, Francesca Olsen, and Neil Johnson. Efficiency of energy 

transfer in a light-harvesting system under quantum coherence. Phys. Rev. B, 78:085115, 

2008. 

[10] Masoud Mohseni, Patrick Rebentrost, Seth Lloyd, and Alan Aspuru-Guzik. Environment- 
assisted quantum walks in photosynthetic energy transfer. J. Chem. Phys., 129:174106, 2008. 

[11] Seth Lloyd. The quantum Goldilocks effect: on the convergence of timescales in quantum 
transport. | arXiv: 1111. 4 982\ 2011. 

[12] Jianshu Cao and Robert J Silbey. Optimization of Exciton Trapping in Energy Transfer 
Processes. J. Phys. Chem. A, 113:13825, 2009. 

[13] Alejandro Perdomo, Leslie Vogt, Ali Najmaie, and Alan Aspuru-Guzik. Engineering directed 
excitonic energy transfer. Appl. Phys. Lett., 96:093114, 2010. 

[14] S M Vlaming, V A Malyshev, and J Knoester. Nonmonotonic energy harvesting efficiency in 
biased exciton chains. J. Chem. Phys., 127:154719, 2007. 

[15] Masoud Mohseni, Alireza Shabani, Seth Lloyd, and Herschel Rabitz. Optimal and robust energy 
transport in light-harvesting complexes: (II) A quantum interplay of multichromophoric 
geometries and environmental interactions. arXiv: 1 104-4812 2011. 

[16] Lorenzo Campos Venuti and Paolo Zanardi. Excitation transfer through open quantum 
networks: Three basic mechanisms. Phys. Rev. B, 84:134206, 2011. 

[17] Stephan Hoyer, Mohan Sarovar, and K Birgitta Whaley. Limits of quantum speedup in 
photosynthetic light harvesting. New J. Phys., 12:065041, 2010. 



Environment- assisted quantum transport in ordered systems 12 

[18] P Nalbach, J Eckel and M Thorwart. Quantum coherent biomolecular energy transfer with 

spatially correlated fluctuations. New J. Phys., 12:065043, 2010. 
[19] P Nalbach, D Braun and M Thorwart. Exciton transfer dynamics and quantumness of energy 

transfer in the Fenna-Matthews-Olson complex. Phys. Rev. E, 84:041926, 2011. 
[20] Torsten Scholak, Fernando de Melo, Thomas Wellens, Florian Mintert, and Andreas Buchleitner. 

Efficient and coherent excitation transfer across disordered molecular networks. Phys. Rev. 

E, 83:021912, 2011. 

[21] Torsten Scholak, Thomas Wellens and Andreas Buchleitner. Optimal networks for excitonic 
energy transport. J. Phys. B: At. Mol. Opt. Phys., 44:184012, 2011. 

[22] M Esposito and P Gaspard. Exactly Solvable Model of Quantum Diffusion. J. Stat. Phys., 
121:463, 2005. 

[23] M Esposito and P Gaspard. Emergence of diffusion in finite quantum systems. Phys. Rev. B, 
71:214302, 2005. 

[24] Marko Znidaric. Exact solution for a diffusive nonequilibrium steady state of an open quantum 

chain. J. Stat. Mech., page L05002, 2010. 
[25] Marko Znidaric. Solvable quantum nonequilibrium model exhibiting a phase transition and a 

matrix product representation. Phys. Rev. E, 83:011108, 2011. 
[26] E Schwarzer and H Haken. The moments of the coupled coherent and incoherent motion of 

excitons. Phys. Lett. A, 42:317, 1972. 
[27] Viv Kendon. Decoherence in quantum walks — a review. Math. Struct. Comp. Set., 17:1, 2007. 



